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Abstract 

We present examples and diagrams illustrating the proofs appearing 
in [6j, to which this paper is meant to be an appendix. We show how 
matrix calculations may be represented by topological surface gluings, 
which may in turn be represented as permutations. 

1 Introduction 

Random matrix calculations can often be represented by topological con- 
structions (see [I], 5 , 9 1 ) . These topological constructions can be represented 
by permutations, as in [IJ[5j[8], which are often more convenient for the 
proofs. This document is intended as an appendix to [6] (which we will refer 
to as the main paper), in which we attempt to provide some intuition for 
how the topological constructions are represented in the permutations, and 
describe the pictures which motivated the various proofs. We illustrate with 
a number of examples. 

In Section [2j we show how a matrix calculation involving Gaussian ran- 
dom matrices can be interpreted in terms of face gluings. As pointed out 
in the exercises to Chapter 3, calculations with real matrices include 
nonorientable surfaces. In Section [3] we provide some intuition for how the 
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Figure 1: The faces used to calculate the expected value of a product of 
traces of products of Gaussian random matrices given in the example. 

permutations, and in particular premaps, may be used to represent these face 
gluings, even in the nonorientable case, and thus used in the matrix calcu- 
lations. In Section [4] we suggest pictures for the Ginibre and GOE matrices 
(Wishart matrices are discussed in the main paper). We suggest intuition 
for further quantities calculated for the matrices: in Section [5] we describe 
how cumulants are given by connected surfaces, in Section [6] how expressions 
with centred terms are given by diagrams where the regions corresponding 
to those terms are not disconnected and in Section [7] we descibe asymptotic 
quantities in terms of topology and noncrossing diagrams. In Section [8j we 
describe how the conditions discussed in Section [6] and Section [7] constrain 
the possible diagrams contributing to the cumulants of traces of products of 
centred matrices, those considered in first- and second-order freeness. 

2 Matrix Calculations as Face Gluings 

As an example of how a topological surface may be used in a matrix cal- 
culation, let fij be independent N (0, 1) random variables, and let X be an 
M x N matrix with Xij = ^r^fij- Let Y\, . . . , Y n be random matrices inde- 
pendent from X, and let 7 G S n and e : [n] — > {1, —1}- In Lemma 3.4 of the 
main paper, we calculate expressions of the form 

E(tr 7 (x^Y 1 ,...,X^Y n )). 

We outline the topological constructions corresponding to the calculations 
in the proof. 

Lemma 3.1 of the main paper suggests that the constraints on the indices 
appearing in a trace have a cyclic structure. We construct a face for each 
trace in the expression, with the edges representing the terms from X and 
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Figure 2: An untwisted edge identification occurring when the term f L + L - 

from the fifth random matrix term X T is constrained to be equal to the 
term / + - from the eighth matrix term X. 

b 8 l 8 




Figure 3: A twisted edge identification occurring when the term f + - from 
the first random matrix term X is constrained to be equal to the term f + - 
from the seventh matrix term X. 

X T ', and the vertices the cyclically in the order they appear in the cycles 
of 7. In Figure [T] we show the faces we would construct to calculate 

E (tr (XY 1 XY 2 X T Y 3 Y 4 X T Y 5 ) tr (X T Y 6 XY 7 XY 8 )) . 

(We choose indices so that the fcth occurrence of the matrix X appears with 
indices it and CZ , which means the order they appear in depends on whether 
this occurrence of X appears with a transpose or not. The arrows designate 
this direction. The shared indices are shown between the matrices sharing 
them.) We apply the Wick formula (see [4j, p. 164) to the random variables. 
The entries of X are independent random variables, so for a given pairing, 
the indices on the paired random matrix entries must be equal, which results 
in another set of constraints on the indices. We represent these constraints 
by gluing the paired edges together so that the indices which must be equal 
are next to each other. In this way we form a surface. 

We note that if one term is from an X and another from X T (as in 
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Figure 4: An example of a gluing of the faces shown in Figure [TJ 




Figure 5: A vertex which appears in the surface gluing shown in Figure |4j 
shown from both sides. 

Figure [2]) , the ordering of the indices means that the edge identification is 
untwisted. However, if both are from an X or an X T (as in Fi gure o ) , 
then the edge identification is twisted. This is where differences from the 
complex case begin to appear. If Z is a complex Gaussian random variable, 
then E (Z 2 ) = E (z 2 ) = 0, while E (ZZ) = 1. Thus one term of a paired 

term must be from an X and the other from an X* , so the twisted edge 
identifications do not happen. This excludes nonorientable surfaces, as well 
as the need to consider other relative orientations of the faces in diagrams 
drawn in the plane. It is the latter that leads to the differences between the 
asymptotic fluctuations of the real and complex matrices. 

We show a surface gluing in Figure |4j The resulting surface will have 
a number of vertices. The constraints on the indices around each are those 
that would appear in a trace, so each vertex corresponds to a trace of the 
Yk matrices (see Figure [5]) . We note that some of the corners of the faces 
may be upside down relative to the others. Since this reverses the order of 
the indices as they appear, we replace that with Yu . (In general, we can 
imagine that the transposes of the matrices are written on the backs of the 
faces.) We note that if we consider the same vertex from the opposite side, 
the matrix product is the transpose of that on the original side, so the trace 
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Figure 6: A simple map, as described in Example |3.1[ 



is unchanged. In this case, the vertex contributes 

Tr (FxlfFelflf) = Tr (Y^Y^Y?) . 

If we are considering an expression in normalized traces, a more natural 
quantity when we consider large N limits, then each trace contributes a 
factor of N. Thus, the highest order terms are those with the most vertices, 
that is, those with the largest Euler characteristics. 



3 Cartography on Unoriented Surfaces 

The face gluings (and hence the constraints on the indices) may be expressed 
in terms of permutations. However, the encoding used in [IppC) depends on 



the surface having an orientation. We outline how we may similarly encode 
an unoriented face gluing. 

If we have an oriented surface gluing, the ends of the edges have a cyclic 
ordering around the vertex they appear at. We can define a permutation 
on the edge-ends a (for sommet, French for vertex) with a cycle for each 
vertex in which the edge-ends at that vertex appear in their cyclic order, 
counterclockwise. Similarly, we can define a permutation a (for arete, edge) 
whose cycles contain the two ends of each edge. Then the cycles of the 
permutation ip := cr _1 Q! _1 enumerate the edges appearing around each face 
of the surface gluing (with the convention that the edge is listed by its first 
edge-end encountered in a counterclockwise direction around the face; the 
other face sharing that edge will list the other edge-end) . A graph with this 
sort of surface embedding information, or this description of it in terms of 
permutations, is often called a map. 

Example 3.1. Consider the map in Figure [6} The edge-ends are labelled by 
the integers in [6]. Then a = (1, 2, 3, 4) (5, 6), and a = (1, 5) (2, 3) (4, 6). 
We calculate that g~ x cT = (1,6,3) (2) (4,5). The first cycle corresponds 
to the outside face, the second to the loop on the left, and the third to the 
two-edged face on the right. We note that the face with two edges lists its 
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Figure 7: A hypermap map, as described in Example 



3.2 



edges as 4 and 5, while these edges are referred to as 6 and 1, respectively, 
when they are listed by the outside edge. 

It is also possible to represent hypermaps in this way. A hypermap is a 
map which may have hyperedges, that is, edges which may have any positive 
integer number of ends, rather than just two. The ends of a hyperedge 
must also have a cyclic ordering representing their embedding in the surface. 
Hyperedges may be visualized as either a separate set of vertices or faces, 
with the edge-ends becoming the edges. 

Example 3.2. In Figure [7| we show a hypermap. The hyperedges are shown 
as outlined shapes connecting them to the black circles representing the 
vertices. The edge-ends are labelled by integers in [7]. We have 

a = (1,2, 3) (4,5) (6,7) 

and 

a = (1,6, 5) (2, 7, 3) (4). 

We can calculate that 

a^oT 1 = (1,4,5,7) (2) (3,6). 

The first cycle corresponds to the outside face, the second to the region on 
the left with one vertex, and the third to the region with two vertices. 

While both examples are planar in order to have clearer diagrams, it is 
possible to use these descriptions for maps on higher genus surfaces. How- 
ever, the description depends on the surface having a consistently defined 
orientation. In order to describe the nonorientable surfaces which may ap- 
pear in the real matrix problems, we extend this construction to the nonori- 
entable case. 

Any surface has an orientable covering space, which we may define by 
constructing two preimages of every point corresponding to the two pos- 
sible orientations. Points are close in the covering space if they are close 
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Figure 8: The orientable covering space of the gluing shown in Figure |4j 

in the original space with the same orientation. For more details on this 
construction, see j2], pages 234-235. We can think of the covering space as 
the surface experienced by someone on the surface, rather than within it. 

For each face of our original surface, we construct two preimage faces, 
which we can think of as the front and the back. We label the latter with 
negative signs. When we glue edges in an untwisted way, we glue the cor- 
responding edges on the fronts and the corresponding edges on the backs, 
while when we glue edges in a twisted way, we glue the corresponding edge 
of the front to the corresponding edge of the back, and vice versa. If the 
original surface was nonorientable, it is possible to travel from the front to 
the corresponding point on the back within the covering space. If the origi- 
nal surface was orientable, we obtain two disconnected copies of the original 
surface. In Figure [8j we show the cover constructed in this way from the 
gluing shown in Figure |4j 

In the diagrams representing the matrix calculations in this work, we will 
generally begin with faces and edge or hyperedge information, and calculate 
vertex information, but this is simply the dual problem. If we begin with 
face information described by a permutation 7, we construct two preimage 
faces for each face: the fronts described by permutation 7+, and those on the 
back described by 7Z 1 , since positive orientation is in the opposite direction 
on the back. We note that 7+7Z 1 , which represents all faces of the covering 
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space, satisfies the definition for a premap. In the example shown in Figure[4j 
we have 

7+ = 7 = (1,2,3,4,5) (6,7,8), 
7 _ = (-1, -2, -3, -4, -5) (-6, -7, -8) , 

and 

7+7I 1 = (1, 2, 3, 4, 5) (6, 7, 8) (-5, -4, -3, -2, -1) (-8, -7, -6) . 

The hyperedges may connect front and back, and therefore a cycle may 
contain both positive and negative integers. However, for any cycle, we 
should be able to find the same cycle from the other side of the surface: a 
distinct cycle which contains the integers with sign reversed and in opposite 
order; thus, the permutation ir representing the hyperedges should also be 
a premap. See [8] for a similar construction. 

We determine the vertex information in a similar manner. Since the 
matrix Y/% appears counterclockwise of the fcth occurrence of X or X T on 
the front but clockwise of it on the back, the convention for which edge-end 
corresponds to which side of an edge means that the permutation 7_7^ 1 7r~ 1 
will list 7_ (—k) instead of —k for any negative integer. We find that we must 
conjugate the permutation by 7Z , giving us the permutation 7^ 1 7r~ 1 7_. 
The traces read off of the vertices are read clockwise, so we will be using the 
inverse of the permutation describing the vertices: jZ 1 K r y + . 

Example 3.3. Let 7 = (1, 2, 3, 4, 5) (6, 7) describe the faces of a nonorientable 
map, and let vr = (1, 7, -6) (6, -7, -1) (2, 5, -3) (3, -5, -2) (4) (-4). We 
show this construction in Figure [9j and represent the change in sign in a 
cycle by a twisting of the hyperedge (this represents that the edge is glued 
in reverse direction). We calculate that 

7+ 1 tt- 1 7_ = (1, -6, 7, 5) (-5, -7, 6, -1) (2, -3, -4) (4, 3, -2) . 

In the matrix calculations, we would use the inverse of this permutation. 
We may trace the vertices in the diagram by following the boundaries of the 
faces and hyperedges, respecting twists. The vertex (1, —6, 7, 5) is shown in 
the diagram. 

We may think of the hyperedges as either faces or vertices and the edge- 
ends as conventional edges. In Figure [7J hyperedges (in grey) may be seen 
as faces alternating with the original faces (in white), while the boundaries 
of the hyperedges (in black) may be seen as the edges (with the convention 
that each edge-end corresponds to the boundary clockwise from it relative to 
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Figure 9: An example of a nonorientable hypermap, discussed in Exam- 
ple 3.3 
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Figure 10: The hypermap in Figure [7] where the hyperedges are seen as a 
set of vertices, shown in white, which alternate with the original vertices, 
and the edge-ends are seen as the edges. 



9 



the hyperedge). Figure 10 shows the graph when we think of the hyperedges 
as vertices (white circles) which alternate with the original vertices (black 
circles), and again the edge-ends (black lines) are the edges. In either case, 
we get the following formula for the Euler characteristic (Definition 3.6 in 
the main paper): 

X (7, tt) := # (7+7I 1 ) /2 + # (tt) /2 + # frl 1 *"^) /2 - |/| . 

Example 3.4. We return to the example shown in Figure [4} We have 

vr = (1, -7) (-1, 7) (2, -4) (-2, 4) (3, -6) (-3, 6) (5, 8) (-5, -8) . 

If the edge identification is untwisted, then in the covering space the front is 
glued to the front and back to the back, so the integers with the same sign 
appear in the same cycle. If the edge identification is twisted, then the front 
is glued to the back and the back to the front, so integers with opposite 
signs appear in the same cycle. We calculate that 

71^7+ = (1, -3, 6, -5, -7) (7, 5, -6, 3, -1) (2, -8, -4) (4, 8, -2) 

and 

7r 1 7T7+/2 = (1, -3, 6, -5, -7) (2, -8, -4) . 

The vertex represented by the first cycle is the one shown in Figure [5} 

We calculate that x ( 7 > 7r ) = 2 + 4 + 2 — 8 = 0. Since the surface is 
connected and nonorientable, it must be a Klein bottle. 

4 The Matrix Models 

The expansion satisfied by all the matrix models in the main paper: 

E(tr 7 (4f »Yi,.-., X« n)) Y n 

= £ iV x(7,4^)-2#( 7 ) /c (?r) E (tr 7ll ^_ l47+/2 (Yi, . . . , Y n 
nePM c (±[n]) 

may be interpreted as a set of rules for constructing surfaces out of the given 
faces. 

In the case of Ginibre matrices, we may interpret the expression as a 
sum over all surface gluings which are consistent with directions built into 
the edges of the faces (counterclockwise if the matrix appears untransposed, 
and clockwise if it appears transposed). We give an example: 
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Figure 11: The faces constructed in order to compute the expression in 
Example 4.1 with a possible edge identification marked. 



Example 4.1. Let Z be a real Ginibre matrix. If we wish to calculate the 
expression 

E (tr (ZY 1 ZY 2 Z T Y 3 Z T Y 4 ) tr (ZY 5 Z T Y 6 ZY 7 Z T Y 8 )) , 



we construct the faces shown in Figure 11 We have 7 = (1, 2, 3, 4) (5, 6, 7, 8) 
and S e = (1) (-1) (2) (-2) (3, -3) (4, -4) (5) (-5) (6, -6) (7) (-7) (8, -8). 
The pairing shown is 

Tr = 5 eP 5p5 e = (1, -7) (-1, 7) (2, 3) (-2, -3) (4, -6) (-4, 6) (5, 8) (-5, -8) 

for p = (1, 7) (2, 3) (4, 6) (5, 8). Paired integers of the same sign represent an 
untwisted edge identification and paired integers of opposite signs represent 
a twisted one. We compute that 

71^7+ = (1, 3, -5, -7) (7, 5, -3, -1) (2) (-2) (4, -6) (6, -4) (8) (-8) 

and that xil^) = 2 + 4 + 4 — 8 = 2, so the example is a sphere (since it is 
connected), and contributes the term 

E (tr (Y 1 y 3 n T *7 T ) tr (Y 2 ) tr (Y 4 Y 6 T ) tr (Y 8 )) N~ 2 . 

The expression for GOE matrices may be interpreted as a sum over all 
edge gluings, in either direction. 

Example 4.2. If we wish to calculate E (tr (r 2 )), where T is a GOE matrix, 
we construct a face with two edges. We may identify the two edges in an 
untwisted way, making the surface a sphere (Euler characteristic 2), or in 
a twisted way, making the surface a projective plane (Euler characteristic 
1) (see Figure 12). These two surfaces correspond to the two elements of 
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Figure 12: The two surfaces corresponding to the two terms contributing to 
the expected value calculated in Example |4.2| 




Figure 13: A surface with two connected components. 

PM(±[2])nP 2 (±[2]), that is, (1,2) (-1,-2) and (1,-2) (-1,2) respec- 
tively. We calculate that 

E (tr (T 2 )) = l + JV -1 . 

If we wish to calculate E (tr (TY1T2Y2)), where Y\ and Y2 are random 
matrices independent from T, we get 

E (tr (TY X T 2 Y 2 )) = E (Y x ) (Y 2 ) + E (Y 1 Y 2 T ) N' 1 ; 

the matrices Y\ and Y 2 appear in different vertices in the sphere, but appear 
in the same vertex but opposite orientation in the projective plane. 

An example calculation involving Wishart matrices and another involv- 
ing several independent matrices are given in the main paper. 

5 Cumulants 

Cumulants are sums over only connected surfaces: the moment-cumulant 
formula organizes the terms by which faces are part of a connected compo- 
nent. The terms which appear in the term corresponding to partition (on 
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elements corresonding to the traces, that is, the faces) tt are those where the 
faces each block of tt form a connected component. The surface constructed 



in Figure 13 for the calculation of 

E (tr (zW^) tr {\V[ X3) Z T Z T ) tr (w^ 6) Z T W^ 8 Vf 9 ^ 
would appear in the term 

k 2 (tr (zwf 2) ) ,tr (wf s) Z T wf ^I^)) h (tr (w[ M) Z T Z T ' 
These terms account for all the disconnected terms, so only the terms con- 



necting all three faces, such as the one in Figure 14, contribute to the cu- 
mulant of the three traces: 

h (tr (ZW^ 2) ) ,tr (w[ H) Z T Z T ) ,tr (w^ X(i) Z T W^ Xs Vf 9 ^ 

6 Centred Terms and the Principle of Inclusion 
and Exclusion 

When we expand a cumulant of traces of products of centred terms, 
k r (tr (ii • • • A P A , . . . , tr (A Pl+ ... +Pr _ 1+ i ■■■A p 

= E (-i) |iq II E ( tr (^)) 

KC[p] k&K 

k r tr ]J A k ,...,tr [] A k 

\ \k&[l, Pl }\K J \ke[ Pl +-+Pr-l + l,p]\K 

we can think of the terms with indices in subset K, which appear in the first 
product, as being required to have a diagram not connected to any other 
(that is, one that would contribute to E (tr (Ak)). Other terms may also be 
disconnected. 

For example, consider the following cumulant of traces of products of 
centred terms: 



{tr [(Z - E (tr (Z))) {w^ M) - E (tr {w^ M) 



tr 
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Figure 14: A term with some connected terms and some disconnected terms. 



(Note that we have decided to consider the two Ginibre matrices in the sec- 



ond trace as a single term, so we have centred their product.) In Figure 14 



the second term in the first trace W^ 2 \ the second term in the second trace 
Z T Z T , and the first and third terms in the third trace and W2 are all 
disconnected. (Again, even though the two Ginibre matrices in the second 
trace are connected to each other, this term is not connected to any other 
term, so it is considered disconnected.) The set of indices of disconnected 
terms is then {2,4,5,7} C [8]). It then appears with sign (—1)'^' in g{K) 
for any K with K C {2, 4, 5, 7}, but only appears in / ({2, 4, 5, 7}). 

We can think of the alternating sign (—1)'^' as subtracting out all 
terms with one disconnected term (with over-counting), then compensat- 
ing by adding back those with two disconnected terms (again with over- 
counting), then again compensating by subtracting those with three discon- 
nected terms, and so on. We are left with only those terms where no term 
is disconnected. 



7 Genus Expansions and Asymptotics 

Moments and cumulants of each of the matrix models we consider can be put 
in a form similar to what is called a genus expansion, or more accurately 
in this case an Euler characteristic expansion. Such a form is a sum of 
terms corresponding to surfaces, where the order of a term in the size of the 
matrix N depends only on the Euler characteristic of the surface. Highest 
order terms correspond to highest Euler characteristic, either spheres or 
collections of spheres. Gluing diagrams for spheres or collections of spheres 
can also be drawn as noncrossing diagrams, since crossings require handles 
to accommodate them. In the real case, in which there are often two possible 
directions in which edges may be identified, the edge identifications must not 
be twisted, since the nonorientable surfaces have lower Euler characterstic 
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Figure 15: An example of a sphere constructed from the faces shown in 
Figure [T] The face on the right has been flipped over, so it is the back that 
is visible. 



than the sphere. However, we must now consider all relative orientations of 
the faces or circles. Figure [15] shows an example of a noncrossing diagram 
on the faces from Figure [T] which results in a sphere. However, the triangle 
is flipped over (that is, its orientation relative to that of the pentagon is 
reversed) in order to demonstrate that the diagram is noncrossing. It is in 
this orientation that the edge-identifications are untwisted. 

The classification theorem for compact connected surfaces states that 
all connected, compact surfaces are either spheres (Euler characteristic 2), 
connected sums of n tori (Euler characteristic 2 — 2n), or connected sums 
of n projective planes (Euler characteristic 2 — n) (see, for example, |3j, 
section 21, for a combinatorial proof). The highest order terms of cumulants, 
which correspond to a single sphere, can be drawn as connected noncrossing 
diagrams. In particular, covariances correspond to connected noncrossing 
diagrams on two circles, or annular noncrossing diagrams, which are studied 
in (5l[7l. 



Example 7.1. There are three pairings on 4 points, shown in Figure 16 Two 
are noncrossing (the one on the left and the one on the right). To contribute 
a highest-order term, the edge-identifcations must also be untwisted. If the 
matrices are Ginibre matrices, it must therefore be a *-pairing. For the 4- 



gon on the left in Figure 11 only the centre and right pairings are *-pairings, 



so only one term contributes to the asymptotic moment. We get: 

lim E (tr (ZY 1 ZY 2 Z T Y 3 Z T Y 4 )) = E (tr tr (Y 2 ) tr (Y 4 )) , 
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Figure 16: The three possible pairings on the edges of a square face; the 
ones on the left and right are noncrossing. 



or, if Yk = In for all k, we can calculate by counting diagrams: 

lim E(tr (ZZZ T Z T )) = 1. 



For the 4-gon on the left in Figure 1 1 , the left and right pairings (but not 



the centre one) are *-pairings, so two terms contribute to the asymptotic 
moment. We get: 

lim E(tr(ZY 5 Z T Y 6 ZY 7 Z T Y 8 )) 

= E (tr (Yi) tr (Y 2 Y 4 ) tr (Y 3 )) + E (tr (Y^) tr (Y 2 ) tr (Y 4 )) 

or 

lim Eftr (ZZ T ZZ T )) = 2. 

The asymptotic values of the fluctuations are given by the two-face maps 
appropriate to the matrix model which form a single sphere. These corre- 
spond to the connected noncrossing permutations that can be drawn on 
the two faces. In the real case, we must consider whether a permutation is 
noncrossing on both possible relative orientations of the two faces. 

The fluctuations of the Ginibre matrices are given by the number of 
annular noncrossing *-pairings under both relative orientations of the two 
faces, where the unstarred terms become starred and vice versa when the 



orientation of a face is reversed (compare Figure 17 with Figure 11) 
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Figure 17: The faces and pairing from Figure 11 with the face on the right 
flipped (over a horizontal axis). 



Example 7.2. To calculate the asymptotic covariance of tr 
tr (ZZ T ZZ T ), we construct the faces shown in Figure 
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ZZZ T Z T ) and 
We show the 

same faces with the same pairing in Figure 17 however, we have flipped the 
face on the right over its horizontal axis. On the "back" of the face, the 
transpose of the matrix on the front is shown. We can see that this pairing 
is a *-pairing under this relative orientation, where it is also noncrossing. 

We may calculate the asymptotic covariance by counting the annular 
noncrossing *-pairings on both relative orientations of the faces, which we 
may do by inspection in this case. The pairing must connect the two faces, 
so there must be at least two pairings between the faces, which must pair an 
adjacent pair of edges to another adjacent pair of faces to avoid crossings. 
We find that there are 4 contributing pairings in the relative orientation 



shown in Figure 11 and 4 in the relative orientation shown in Figure 17 
Thus, 

lim k 2 (Tr (ZZZ T Z T ) ,Tr (ZZ T ZZ T )) = 8. 



8 Freeness 

In Lemma 6.1 and Theorem 6.2 of the main paper, we consider expressions 
where the centred terms come from alternating or cyclically alternating al- 
gebras, so the terms cannot be connected to their neighbours. If we wish 
to consider the terms which survive asymptotically, we require that the dia- 
gram be noncrossing. If a term is connected to a non-neighbouring term on 
the same face, this separates the terms between the connected terms from 
all the others. This eventually forces an unconnected term (forbidden by 
the centring of the terms) or a term paired with its neighbour (forbidden 
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Figure 18: An example of a spoke diagram. Dashed lines mark which edges 
belong to a single term, and dotted lines separate the spokes. 
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Figure 19: The spoke on A\ and which corresponds to a noncrossing 
diagram on a disc. 



by the cyclic alternating of the terms). In the one-trace expressions, we 
are left with no possible diagrams, so the quantity vanishes asymptotically, 
meaning that the matrices are asymptotically free. In the two-trace expres- 
sions, a term may instead be connected to a term on the other face. This 
significantly restricts the sorts of diagrams which are possible: specifically, 
the terms must be connected in a sort of spoke arrangement, as in Figure 1 



of the main paper (with an example shown in Figure 18 here). In the real 
case, both relative orientations of the two faces must be considered. 

For a given spoke diagram, the contribution is the sum over all choices 
of all possible diagrams (those which are connected, noncrossing, untwisted, 
and compatible with the matrix model) connecting each pair forming a 
spoke. Since the contribution of each spoke is multiplicative, we can factor 
this as the product over the spokes of the sum over the possible diagrams 
on that spoke. These diagrams will be disc- noncrossing on the two intervals 
comprising the spoke (that is, those that would contribute to the limit value 



of the moment, with an example shown in Figure 19), but we must subtract 
out the disconnected diagrams (the product of a choice of a possible term 
on each interval individually, that is, the product of the moments on each 
interval individually), so we get the contribution 

E (tr (^ (±1) )) -E(tr (A k ))E (tr (^ (±1) )) . 

Example 8.1. Let colour 1 correspond to a Ginibre matrix Z, colour 2 cor- 
respond to a GOE matrix T, and colour 3 correspond to a Wishart matrix 
W := X T X, X as in the definition. Let A x = ZZZ T , A 2 = T 2 , A 3 = W 2 , 
B\ = W 5 , B2 = ZZZ 1 ', and B% = T 4 . If we consider the reversed spoke di- 
agram with k = 1 (corresponding to the bottom-centre diagram in Figure 1 



of the main paper), Figure 18 shows a possible diagram 
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By inspection, we can see that this is the only spoke diagram that will 
accommodate a diagram on A\ and B%, the two members of the algebra 
generated by Z. Again by inspection, we can see that there is only one 
other diagram on the face shown in Figure 19, and it also connects A\ to 
B\, so the contribution of this spoke is 

E (tr (ZZZ T (ZZZ T ) T )) - E (tr (ZZZ T )) E (tr ((ZZZ T ) T )) = 2. 



Using the results cited in Remark 5.13 of the main paper, we may calculate 
that: 

E (tr (T 6 )) - E (tr (T 2 )) E (tr (T 4 )) =5-1-2 = 3 

and 

E (tr (W 7 )) - E (tr (W 2 )) E (tr (W 5 )) = 429 - 2 • 42 = 345, 

so 

Jiml! 2 (Tr (A 1 A 2 A 3 S \ ,Tr (b^B^ = 2070. 
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